Mon. Not. R. Astron. Soc. 000, [TH5] (2013) Printed 22 February 2013 (MN M£X style file v2.2) 



On the dynamical stability of the proposed planetary 
system orbiting NSVS 14256825 



cn 
o 

(N 
X> 

IX, 
(N 

6 

o3 



> 

m 

O 

m 



Robert A. Wittenmyer 1 *, J. Horner 1 and J. P. Marshall 2 

1 Department of Astrophysics and Optics, School of Physics, University of New South Wales, Sydney 2052, Australia 
2 Departmento Fisica Teorica, Universidad Autonoma de Madrid, Cantoblanco, 28049, Madrid, Spain 



Accepted Received ; in original form 



ABSTRACT 

We present a detailed dynamical analysis of the orbital stability of the two 
circumbinary planets recently proposed to orbit the evolved eclipsing binary star 
system NSVS 14256825. As is the case for other recently proposed circumbinary 
planetary systems detected through the timing of mutual eclipses between the central 
binary stars, the proposed planets do not stand up to dynamical scrutiny. The 
proposed orbits for the two planets are extremely unstable on timescales of less than 
a thousand years, regardless of the mutual inclination between the planetary orbits. 

For the scenario where the planetary orbits are coplanar, a small region of 
moderate stability was observed, featuring orbits that were somewhat protected from 
destabilisation by the influence of mutual 2:1 mean- motion resonance between the 
orbits of the planets. Even in this stable region, however, the systems tested typically 
only survived on timescales of order 1 million years, far shorter than the age of the 
system. 

Our results suggest that, if there are planets in the NSVS 14256825 system, they 
must move on orbits dramatically different to those proposed in the discovery work. 
More observations are clearly critically required in order to constrain the nature of 
the suggested orbital bodies. 
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1 INTRODUCTION 

The intriguing possibility of circumbinary planets has been 
a feature of science-fiction for decades. Theoretical studies 
of planetary f ormation have shown that such systems are in - 
deed possible (|Quintana fc Lissauer 2006; G ong et al.ll2013h . 
and that they can be dyna mically stable (jKubala et al .11 19931 : 
iHolman fc Wiegertl Il999h . but their secure detection re- 
mained quite elusive, even in this present age of accelerating 
exoplanetary discovery. Finally, in late 2011, the first unam- 
biguous example of a circum binary planet was announced 
by the Kepler science team (|Dovle et alj l201lh . This dis- 
covery was quickly followed by the first multiple- planet 
circumbinary system, Kepler-47 JOrosz et al. 2012alk and 
two additional c ircumbinary systems (|Orosz et al. 2012b: 
IWelsh et al.ll2012h . These exciting results have yielded a first 
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estimate that ~1% of close binary stars host such circumbi- 
nary planets. 

In addition to the Kepler discoveries, there is a growing 
body of literature claiming the detection of circumbinary 
planets orbiting post-common-envelope host stars. Such ob- 
jects can, in theory, be detected by light-travel-time vari- 
ations: as the postulated planetary companions orbit the 
central stars, they cause those stars to move back and forth 
as they orbit around the system's centre of mass. As a re- 
sult, the distance between the Earth and the host stars 
varies as a function of time, meaning that the light from 
the stars must sometimes travel further to reach us than at 
other times. This effect results in measurable variations in 
the timing of mutual eclipse events between the two stars 
that can be measured from the Earth. The stars proposed 
as hosts of circumbinary planets detected in this manner in- 
clude cataclysmic variables (e.g. UZFor, Potter et al. 2011), 
pre-cataclysmic variables (NNSer, Beuermann et al. 2010), 
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and detached subdwarf-M dwarf (e.g. HWVir, Beuermann 
et al. 2012) or subd warf-white dwarf (e.g. RRCae, Qian 
et al. 2012) binaries. IZorotovic fc Schreiberl (120131 ) provide 
a complete summary of the 12 proposed post-common- 
envelope planetary systems in their Table 2. 

While the profusion of circumbinary planets claimed 
to orbit post-common-envelope binaries would suggest that 
such systems are extremely common, all is not as it seems. 
When subjected to rigorous dynamical testing, a great many 
of these proposed planetary systems have been found to be 
simply unfeasible. That is, the orbits of the proposed plan- 
ets are such that, if they were truly planets, they would 
rapidly experience significant and ultimately catastrophic 
mutual interactions (ejections or collisions). Recently, sev- 
eral of the proposed circumbinary multiple-planet systems 
have been tested in this manner. Highly detailed dynam- 
ical maps, showing the lifetimes of a wide range of or- 
bital configurations for the planets in question, demon- 
strate that a number of these candidate systems are dy- 
namically unstable on timescales of J$ 10 4 yr. Some pro- 
posed circumb inary systems investigated in this way in- 
clude HU Apr (IHorner et al.l[201ll; I Wittenmyer et al.ll2012al ; 



Table 1. Data on the NSVS 1425 system from lAlmeida etaU 
ll2012bh . 



Iffinse et alll2012h. NN Ser ^Horner et alj|2012bl ), and HW 
Vir ( Horner et al.ll20l"2T ) . The same dynamical mapping tech- 
niques have also been applied to multiple-planet systems 
discovered by radial velocity studies, in order to check their 
long-term stability l|Wittenmver et al.ll2012bl) and to asse ss 
the role of low-order resonances (|Robertson et al.l l2012aH bl; 
IWittenmver et all [20123 ). 

NSVS 14256825 (hereafter NSVS 1425) is an eclips- 
ing binary consisting of a subdwarf OB star and an M 
dwarf orbiting one ano ther with a period of 0.110374 days 
l|Almeida et al.1 l2012al ). The combination of a hot subd- 
warf and a late-type dwarf produces significant reflection 
effects in the light curves, and such systems are known 
as "H W Vir" system s , of w hich only 10 are currently 
known (|Almeida et alj l2012bh . Sinc e its discovery in the 
Northern Sky Variability Survey (IWozniak et 
eclipse timings have bee n repo rted by IWils et 
iKilkennv fc Koenl (|2012l ), and iBeuermann et al.l (l2012al). 



A lin ear p eriod change was noted by Kilkenny fc Koenl 
i|2012l) . and B euermann et al.l (|2012al ) reported a cyclic pe- 
riod change, which they attributed to the presence of a 
~ 12 Mjup pla net with an un c onstrai ned period P <; 20 yr. 
Most recently, lAlmeida et al.l l|2012bl ) presented additional 
eclipse timings and, by combining all available timings, fit- 
ted two periodicities which were then attributed to the in- 
fluence of two unseen circumbinary planets. The reported 
parameters for the candidate planets are given in Table [T] 
In this work, we bring our well-tested dynamical mapping 
techniques (Section 2) to bear on the candidate NSVS 1425 
planetary system. We determine the dynamical stability of 
the complete ±3<r range of orbital parameters, and we test 
the effect of mutual inclinations between the two planets 
(Section 3). We discuss the results and make our conclu- 
sions in Section 4. 



Parameter 


Inner Planet 


Outer Planet 


Eccentricity 


0.0±0.08 


0.52±0.06 


uj (degrees) 


U.4±7 


97.5±8 


Orbital Period (yrs) 


3.49±0.21 


6.86±0.25 


Orbital Radius (AU) 


1.9±0.3 


2.9±0.6 


Msini (M Jup ) 


2.8±0.3 


8.0±0.8 



withi n the TV-body dynamics package Mercury (|Chambers I 
Il999l) to perform our integrations. We held the initial orbit 
of the in ner planet fixed at its best-fit parameters, as re- 
ported in lAlmeida et all l|2012bl ). and then created 126,075 
test systems. In those test systems, the initial orbit of the 
outer planet was varied systematically in semi-major axis 
a, eccentricity e, periastron argument uj, and mean anomaly 
M, resulting in a 41:r41:rl5:r5 grid of "clones" spaced evenly 
across the 3a range in those parameters (a s given in Ta- 
ble[lj. For the mean anomaly, not reported in lAlmeida et al] 
(|2012br ). we simply tested the entire allowed range 0-360 de- 
grees. We assumed the central stars to be a single point 
mass, an acceptable approximation as the binary separation 
is <C orbital radius of the inner planet, with M = 0.62M© 
l|Beuermann et all I20l2al 1'l We assumed the planets were 
coplanar with each other and had masses equivalent to their 
minimum mass, Msini. We then followed the dynamical 
evolution of each test system for a period of 100 million 
years, and recorded the times at which either of the planets 
was removed from the system. Planets were removed if they 
collided with one another, hit the central body, or reached 
a barycentric distance of 10 AU. 

In addition to this highly detailed "nominal" run, we 
examined the effects of mutual inclinations between the two 
planets by ru nning further suites of simulations at lower res - 
olution (e.g. IHorner et all 1201 ll ; IWittenmver et al.ll2012d ). 
We varied the initial orbit of the outer planet as above, re- 
sulting in a 21a;21a;5a;5 grid in a, e, uj, and M, respectively 
(for a total of 11,025 test systems). We ran five such scenar- 
ios, with the inclination between the two planets set at 5, 15, 
45, 135, and 180 degrees (the latter two corresponding to a 
retrograde orbit for the outer planet). Again, the 11,025 test 
systems were allowed to run for 100 million years, or until 
the system was destabilised due to ejection or collision. 



3 RESULTS 

Figure [T] shows the stability o f the nominal orbits f or the 
NSVS 1425 system as given in lAlmeida et all (|2012bl ). It is 
immedately apparent that the vast majority of the ±3a pa- 
rameter space is extremely unstable, with mean lifetimes of 
less than 1000 yrs. This result is not particularly surprising, 
given that the great majority of solutions tested place the 



2 DYNAMICAL ANALYSIS 

As in our previous wo rk (e.g. iMarshall et al.l l20ld ; 
IHorner et aklhoill . l2012bl ). we used the Hybrid integrator 



1 We note that lAlmeida et al. | d2012al ) recently determined a 
combined mass of 0.528±0.074 for the central stars, 1.2a smaller 
than our assumed value. This small difference would have no ef- 
fect on the interactions between the planets, only on the overall 
scale of the system, and hence does not affect our results. 
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Figure 1. Dynamical stability of the NSVS 1425 system as pro- 
posed by I Almeida et alj (|2012bft . as a function of the semi-major 
axis, a, and eccentricity, e, of the outer planet. The mean lifetime 
of the planetary system (in logio (lifetime/yr)) at a given a — e 
co-ordinate is denoted by the colour of the plot. The lifetime at 
each a — e location is the mean value of 75 separate integrations 
carried out on orbits at that a — e position (testing a combina- 
tion of 15 unique uj values, and 5 unique M values). The nominal 
best-fit orbit for the outer planet is shown as the small red square 
with ± l-cr error bars. 



two planets on mutually crossing orbits - meaning that close 
encounters between the two are a certainty. 

The nominal best-fit orbit of the outer planet does fall in 
a narrow strip of increased stability, featuring mean lifetimes 
~ 10 6 yr. However, as the subdwarf B host star has evolved 
well off the main sequence, one would expect the system to 
be considerably older than this mean lifetim e. While the age s 
of subdwarf B stars are highly uncertain (|Hu et alJ feoiCO. 
even an A-type progenitor would make the system ~0.5 Gyr 
old. Population-synthesis models of the close binary progen- 
itors of these post-common-envelope systems show that the 
initial mass distribution of primary stars that result in subd- 
warf B stars ranges from 1.0-1.8 Mq. The shortest-lived (i.e. 
highest-mass) progenitors would then evolve off the main 
sequence in ~ 2 x 10 9 yr. Hence, any planets which ex- 
isted before the common-envelope stage would be expected 
to demonstrate dynamical stability on a timescale far longer 
than that exhibited by the NSVS 1425 system. 

The observed strip of moderate stability is attributed 
to the protective influence of the mutual 2:1 mean- motion 
resonance between the two proposed planets. Such resonant 
protection is well known, both in our own Sola r system (e.g. 
iHorner fe Lvka wka 2010; Lykawka et aL 2011 ), and in exo- 



planet ary science (e.g. iRivera et aLlbOloT l Robertson et al.l 
l2012al lbTl. Indeed, it is interesting to compare our results with 
those obtai ned for the proposed planets around HU Aqr (e.g. 
Figure 1 of IHorner et al.ll20TlT ). In both cases, the proposed 
outer planet lies on a highly eccentric best-fit orbit (e ~ 0.5) 
that crosses that of the innermost planet and the two planets 
are close to mutual 2:1 mean-motion resonance. Addition- 
ally, the proposed system proves to be dynamically unstable 
on typical timescales of hundreds or thousands of years. 

We then considered the possibility of non-zero mu- 
tual inclinations between the two candidate planets in the 
NSVS 1425 system. For a significantly interacting planetary 
system, retrograde orbits can provide a dynamically sta- 



ble configuration (e.g. [ Eberle fc Cunt j 2010l: IHorner et all 
l201ll : lQuarles et al.ll2012l : lMorais fc Giuppone I2012T ). Toex- 
plore the dynamical stability of systems resulting from such 
non-coplanar scenarios, we performed a second suite of sim- 
ulations at lower resolution, as described in Section 2. The 
results of the five mutually-inclined scenarios are shown in 
Figure [2] along with the nominal coplanar scenario from 
Figure [1] We see that mutually-inclined scenarios are even 
more unstable than the coplanar scenario. Even for the 180- 
degree configuration (i.e. retrograde and coplanar, panel f), 
the only region of stability appears in the lower right, at 
the lowest eccentricities and largest semi-maj or axis for the 
outer planet. As was the case for HU Aquarii (|Horner et al.l 
2011), the highly stable region is restricted to orbits which 
have periastron distances several Hill radii from the inner 
planet. However, this region lie s well beyond the la u ncer- 
tainties of the orbits derived by I Almeida et all l|2012bl ). 



4 DISCUSSION AND CONCLUSIONS 

W e have shown that the t wo circumbinary planets proposed 
by I Almeida et ail l|2012bl ) are dynamically unfeasible in vir- 
tually any configuration within 3cr of their derived orbital 
parameters. Based on our results, a mechanism other than, 
or in addition to, one (or more) pl anets is needed to ex- 
plain the observed period variations. IZorotovic fc Schreiberl 
(2013) note that about 90% of eclipsing post-common- 
envelope binaries are reported to have timing variations, 
nearly all of which have been attributed to one or more or- 
biting circumbinary planets. Both avenues to such planetary 
systems - survival of the planets through the Asymptotic Gi- 
ant Branch (AGB) and planetary nebula phase, or secondary 
accretion from the ejected envelope of the primary star - are 



possible but highly uncertain (e .g. Postnov fc Prokhorov 



19921: iTavani fc Brookshawlll992l; IPhinnev fc Hansen! Il993l : 



Villaver fc Livid|2007) : [Hansen et alj|2009l; Veras et aljboill : 
Kunitomo et allboittlMustill fc Villaver H2012T ). 



IZorotovic fc Schreiberl ( |2013 ) examined the formation 
and observational statistics of circumbinary planets in order 
to investigate the puzzling trend that fully 90% of post- 
common-envelope eclipsing binaries appear to host plan- 
ets. They note t hat observat i ons o f disk-bearing pre-main- 
sequence stars (|Kraus et al.l [2012) indicate that the cir- 
cumbinary disk lifetime is too short ( < 1 Myr) to form giant 
planets by core accretion. Observational results from Kepler 
also point to a mai n-sequence circum bi nary giant planet 
freque ncy of ~1% l|Welsh et al.1 I2012T ). IWittenmver et all 
|201lf) estimate the frequency of giant planets in 3-6 AU 
orbits to be no more than 37%, based on data on single main- 
sequence stars from the Anglo-Aus tralian Planet Search. 
All of these strands of evidence led IZorotovic fc Schreiberl 
|2013l ) to the conclusion that "virtually all close-compact bi- 
naries. ..are unlikely to be explained by first-generation plan- 
ets." 

Another potential solution to the near-ubiquitous pres- 
ence of planets around evolved binaries is the formation 
of these companions in the post main sequence circumstel- 
lar envelope produced by the subdwarf OB progenitor as 
it evolves through its AGB phase. The amount of material 
cast off by an AGB star (up to 70% of the stellar mass, 
lHabing fc Olofsson lliooi ) is comparable to that of the Min- 
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Figure 2. Dynamical stability for the NSVS 1425 system, for six values of inclination between the two planets. Panels (a) through (f) 
represent mutual inclinations of 0, 5, 15, 45, 135, and 180 degrees, respectively. Panel (a) is a duplicate of Figure [T] reprised here for 
ease of comparison. As in the previous figure, the colour bar represents the log of the mean survival time at each a — e position (testing 
a combination of 5 u values and 5 M values). 



imum Mass Solar Nebula (|Weidenschilling||l977l ) and the 
lifetime of the AGB phase is similar to that of gas-rich pro- 
toplanetary dis cs within which first ge neration planet for- 
mation occurs (Hernand ez et al~l 120071 ). lending plausibil- 
ity to the idea of a second generation of planet formation. 
However, this second generation planet formation scenario 
is highly speculative and requires further scrutiny (and de- 
tailed modelling) before it can be con sidered a viable answer 
l|Akashi fc Soker ll200Sl : |Peretsll2010r i. 

As was the case in earlier studies of propo sed circumbi- 
nary planets in highly evolved systems (e.g. | Horner et al.l 
l201ll ; Iwittenmver et al]|2012al ; iHorner et al.ll2012l ). we find 
that the planetary system proposed in the evolved binary 



system NSVS 14256825 is simply not dynamically feasible. 
Our results suggest that some mechanism other than planets 
must be responsible for the observed eclipse-timing varia- 
tions, and once again highlight the critical importance of 
performing detailed dynamical analyses of potential new 
planetary systems in order to determine whether the pro- 
posed systems make sense. 
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